A model-based insight into the coupling of nitrogen and sulfur cycles in a coastal upwelling system

The biogeochemical cycling in oxygen-minimum zones (OMZs) is dominated by the interactions of microbial nitrogen transformations and, as recently observed in the Chilean upwelling system, also through the energetically less favorable remineralization of sulfate reduction. The latter process is masked, however, by rapid sulfide oxidation, most likely through nitrate reduction. Thus, the cryptic sulfur cycle links with the nitrogen cycle in OMZ settings. Here, we model the physical-chemical water column structure and the observed process rates as driven by formation and sinking of organic detritus, to quantify the nitrogen and sulfur cycles in the Chilean OMZ. A new biogeochemical submodule was developed and coupled to the Regional Ocean Model System (ROMS). The model results generally agree with the observed distribution of reactive species and the measured process rates. Modeled heterotrophic nitrate reduction and sulfate reduction are responsible for 47% and 36%, respectively, of organic remineralization in a 150 m deep zone below mixed layer. Anammox contributes to 61% of the fixed nitrogen lost to N2 gas, while the rest of the loss is through canonical denitrification as a combination of organic matter oxidation by nitrite reduction and sulfide-driven denitrification. Mineralization coupled to heterotrophic nitrate reduction supplies ∼48% of the ammonium required by anammox. Due to active sulfate reduction, model results suggest that sulfide-driven denitrification contributes to 36% of the nitrogen loss as N2 gas. Our model results highlight the importance of considering the coupled nitrogen and sulfur cycle in examining open-ocean anoxic processes under present, past, and future conditions.